%测试圆和椭圆的曲率变差
x=0:0.05:1;
y1=sqrt(1-x.^2);
plot(x,y);
y2=sqrt(1-2.*x.^2);